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ABSTRACT 

We present deep radio continuum observations of the cores identified as deeply embedded 
young stellar objects in the Serpens molecular cloud by the Spitzer c2d programme at a wave- 
. length of 1.8 cm with the Arcminute Microkelvin Imager Large Array (AMI-LA). These ob- 

0^ ' servations have a resolution of w 30arcsec and an average sensitivity of 19 /iJy beam^'. The 

I targets are predominantly Class I sources, and we find the detection rate for Class I objects in 

this sample to be low (18%) compared to that of Class objects (67%), consistent with previ- 
ous works. For detected objects we examine correlations of radio luminosity with bolometric 
luminosity and envelope mass and find that these data support correlations found by previous 
samples, but do not show any indiction of the evolutionary divide hinted at by similar data 
from the Perseus molecular cloud when comparing radio luminosity with envelope mass. We 
conclude that envelope mass provides a better indicator for radio luminosity than bolometric 
' luminosity, based on the distribution of deviations from the two correlations. Combining these 

I new data with archival 3.6 cm flux densities we also examine the spectral indices of these ob- 

jects and find an average spectral index of a] ! = 0.53 ±1.14, consistent with the canonical 
value for a partially optically thick spherical or coUimated stellar wind. However, we caution 
that possible inter-epoch variability limits the usefulness of this value, and such variability is 
supported by our identification of a possible flare in the radio history of Serpens SMM 1 . 
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1 INTRODUCTION 

Recent radio follow-up of the Spitzer catalogue of candidate low- 
luminosity embedded protostars has provided substantial new data 
with which to constrain the correlations observed between radio lu- 
minosity and the other physical parameters of protostellar objects 
in the low-luminosity limit (AMI Consortium: Scaife et al. 2011a; 
201 lb, hereafter Papers I & II). These data have also provided in- 

* We request that any reference to this paper cites "AMI Consortium: Scaife 
etal. 2011". 

t email: ascaife@cp.dias.ie 



formation on a handful of very low-luminosity objects (VeLLOs; 
Young et al. 2004) which are of particular interest because their 
luminosities violate the lower limit set by steady accretion mod- 
els (e.g. Shu 1977), sometimes by over an order of magnitude 
(Evans et al. 2009). This discrepancy provides evidence for alterna- 
tive models of protostellar evolution, such as non-steady accretion 
(Kenyon & Hartmann 1995; Young & Evans 2005; Enoch et al. 
2007). 

This paper will deal predominantly with low mass protostars 
in the earliest stages of their evolution. Class objects represent 
the youngest protostars, generally distinguished from gravitation- 
ally bound pre-stellar cores by the presence of a detectable IR 
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source, a compact radio counterpart, or a collimated molecular out- 
flow (Andre et al. 1999). Such objects have yet to accrete significant 
mass from their surrounding envelopes so that Menv > Me. They 
are very faint at optical and Near Infra-Red wavelengths, but have 
significant sub-millimctro emission. Their SEDs are normally well 
characterized by single blackbodies, typically with 10 7^ ^ 15 K 
(Hatchell et al. 2007), and they possess energetic outflows driven by 
their accretion. Class I objects are considered to be a slightly more 
evolved stage of protostellar evolution where the accretion rate has 
declined, and this is reflected in their less powerful outflows (Bon- 
temps et al. 1996) and a more clearly visible IR counterpart. In 
addition they have developed, or are in the process of developing, 
circumstellar disks. The presence of such disks alters the form of 
the SED, broadening it and flattening the Rayleigh- Jeans tail due to 
dust grain growth. Class II objects (Classical T Tauri stars) repre- 
sent the start of the main sequence phase of stellar evolution, where 
infall is complete and the central object is embedded in an optically 
thick disc. Although protostellar evolution can be broadly encom- 
passed by these definitions there is no strict division between the 
different stages, and classification can be complicated by a number 
of observational factors. 

Unlike the correlations with bolometric and infra-red Ivuninos- 
ity similar to those already known for higher luminosity (> 10^ h;-., ; 
Anglada 1995) protostellar sources derived in Paper I, the observed 
correlation of radio luminosity with envelope mass observed in Pa- 
per n suggested an evolutionary distinction between Class and 
Class I objects. This distinction is potentially explained by a dif- 
ference in the origin of radio emission from protostellar objects 
between these classes, with Class sources producing radio emis- 
sion through an intrinsic mechanism but Class I objects being more 
heavily influenced by local environmental conditions. However the 
poor detection rate for Class I objects at 1.8 cm prevents a conclu- 
sive statement about this difference to be made. 

The Serpens molecular cloud is a highly extinguished, active 
near-by star formation region at a distance of D = 260 ± 10 pc 
(Straizys et al. 1996). The main Serpens cluster (~ 6' extent), some- 
times referred to as "Cluster A" (Harvey et al. 2007), has been ex- 
tensively studied in the infrared and submillimetre (see e.g. Davis 
et al. 1999; Casali et al. 1993; Enoch et al. 2007; Graves et al. 
2010) although there are few surveys of the full extent of the op- 
tical extinction attributed to the larger molecular cloud of lOdeg^ 
(Cambresy et al. 1999). The region of active star formation approxi- 
mately 45' to the south of the Serpens core is known as Serpens/G3- 
G6 (Cohen & Kuhi 1979) or "Cluster B" (Enoch et al. 2007; Har- 
vey et al. 2006). Both the Serpens core and the G3-G6 region have 
extensive multi-frequency observations available, including radio 
surveys at longer wavelengths (Eiroa et al. 2005, Serpens core; 
Djupvik et al. 2006, Serpens/G3-G6). The Serpens cloud contains 
a large proportion of Class I and II young stellar objects (YSOs) 
compared with the Perseus molecular cloud, considered in Paper 
II, although these are found predominantly not in the clusters but 
in the extended cloud (Harvey et al. 2007), supporting the idea that 
the clusters, "A" and "B", are younger than the surrounding cloud. 

In this paper we extend our radio follow-up of the Spitzer cata- 
logue of low-luminosity embedded protostars (Dunham et al. 2008, 
hereafter DCE08) with the Arcminute Microkelvin Imager Large 
Array (AMI-LA; Papers I & II) to cover the Serpens region. In 
§ 2 we describe the sample; in § 3 we briefly introduce the tele- 
scope and describe the data reduction procedures utilised; in § 4 
we describe the results of the observations towards the Cluster A, 
Cluster B and Serpens Filament individually and review the com- 
plementary data at other wavelengths. In § 5 we draw correlations 
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Figure 1. Spitzer 24 fim map of Serpens with the FWHM of the primary 
beam for the AMI-LA fields, listed in Table 3, shown as circles. The loca- 
tions of the regions known as Cluster A and Cluster B are indicated. 



between the characteristic parameters of these objects, and in § 6 
we discuss the implications of these results with reference to de- 
tection statistics for differing protostellar classes, radio spectral in- 
dices and, where measured, outflow momentum fluxes. In § 7 we 
draw our conclusions. 



2 THE SAMPLE 

In this work we observe a sample of candidate deeply embedded 
protostellar cores from the Spitzer catalogue of DCE08. Candidates 
for embedded objects from this catalogue were ranked by DCE08 
as belonging to one of six "Groups", with Group 1 being those 
most likely to be true embedded objects, and Group 6 those least 
likely. The sample of objects observed here comprises all objects 
classified as Groups 1-3 in the Serpens region. Group 1 objects 
are sources which have been confirmed as embedded protostars in- 
dependently of their Spitzer identification. Group 2 objects are not 
confirmed as embedded in a dense core but have ancillary data con- 
firming a molecular outflow, and Group 3 are not confirmed as em- 
bedded and do not have outflow data. We identify the candidates 
from this catalogue by their catalogue number, i.e. DCE08-«;jn, in 
column [1] of Table 1. The coordinates of the individual candidates 
are listed in columns [2] and [3] of Table 1, along with their Group 
and infrared luminosity from Dunham et al. (2008) in colunms [4] 
and [5]. Three sources are found in the Cluster A region (DCE08- 
210, -21 1 c& -215) and four in the Cluster B region (DCE08-208 & 
212-214). The final source, DCE08-202, is outside the cluster re- 
gions but is found within the elongated sub-mm filament identified 
in Enoch et al. (2007) to the north of Cluster B, referred to as the 
Serpens filament. The locations of the sources Usted in Table I are 
shown in Fig. 1. 

The source list in Table 1 has been cross-referenced with the 
larger Spitzer cores to disks catalogue of Evans et al. 2009 (here- 
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inafter EDJ09) in order to provide bolometric luminosities and tem- 
peratures, and 3.6 /im flux densities, as well as the Serpens Spitzer 
catalogue of Harvey et al. (2007; hereafter HMH07). 

Hatchell et al. (2007) used three criteria to distinguish Class I 
from Class 0: firstly that the bolometric temperature, Tboj > 
70 K; secondly that the ratio of bolometric to sub-mm limiinos- 
ity, Lboi/Lsmm > 3000; and thirdly that Ss e^m/Sgso^m > 0.003. 
The final classification for a source was based on the majority re- 
sult from these three criteria. Although no measure of the sub-mm 
luminosity is available for these sources from the literature, all of 
the objects in this sample would be classified as Class I by both of 
the two alternative indicators, giving a majority without the need 
for a third indicator. These indicators have been calculated using 
the physical data from Evans et al. (2009). However Enoch et al. 
(2009) classify three of these objects as Class and it is likely that 
they are borderUne cases as the distinction between Class and 
Class I is not exact. The classifications from Enoch et al. (2009) 
are listed in Column [11] of Table 1 . 

From the correlation of radio luminosity and bolometric lumi- 
nosity derived in Paper II we can predict the expected 16 GHz radio 
flux densities of these sources from the bolometric Ivmiinosities in 
Evans et al. (2009). These predictions are shown in Column [13] of 
Table 1 and have uncertainties of a factor of ~ 2. However, although 
this general trend is observed for detected sources, we note that a 
high proportion of Class I objects are un-detected in the radio. 



3 OBSERVATIONS 

AMI comprises two synthesis arrays, one of ten 3.7 m antennas 
(S A) and one of eight 13m antennas (LA), both sited at the Mullard 
Radio Astronomy Observatory at Lord's Bridge, Cambridge (AMI 
Consortium: Zwart et al. 2008). The telescope observes in the band 
13.5-17.9 GHz with eight 0.75 GHz bandwidth channels. In prac- 
tice, the two lowest frequency channels (1 & 2) are not generally 
used due to a lower response in this frequency range and interfer- 
ence from geostationary satellites. The data in this paper were taken 
with the AMI Large Array (AMI-LA). 

Observations of the 8 objects Hsted in Table 1 were made 
with the AMI-LA between December 2010 and January 2011. 
The eight targets were observed in five separate pointings:DCE08- 
202, DCE08-208 and DCE08-215 were considered sufficiently 
separated to be observed as separate pointings; DCE08-210 and 
DCE08-21 1 are separated by ~ 40" and were therefore observed 
as a single pointing; DCE08-212, 213 & 214 have a maximum sep- 
aration of just over 1 arcmin and were therefore also observed as 
a single pointing. These objects lie within the declination range 
— 2°, outside the standard observing range of the AMI telescope; 
below 5 ~ 15° significant interference from geostationary satel- 
lites is experienced. The effect of violating this declination limit is 
a higher noise level than the telescope specifications due to both the 
heavier data flagging required to remove satellite interference and 
the residual contamination. 

AMI-LA data reduction is performed using the local software 
tool REDUCE. This applies both automatic and manual flags for 
interference, shadowing and hardware errors, Fourier transforms 
the correlator data to synthesize frequency chaimels and performs 
phase and amplitude calibrations before output to disc in uv FITS 
format suitable for imaging in AIPS^ Flux (primary) calibration is 
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Table 2. AMI-LA frequency channels and primary calibrator flux densities 
measured in Jy. 



Channel No. 


3 


4 


5 


6 


7 


8 


Freq. [GHz] 


13.88 


14.63 


15.38 


16.13 


16.88 


17.63 


3C48 


1.85 


1.75 


1.66 


1.58 


1.50 


1.43 


3C286 


3.60 


3.54 


3.42 


3.31 


3.21 


3.11 


3C147 


2.75 


2.62 


2.50 


2.40 


2.30 


2.20 



Table 3. AMI-LA observations. Columns are [1] source number from the 
Dunham et al. (2008) catalogue; [2] AMI-LA flux calibrator; [3] AMI-LA 
phase calibrator; [4] rms noise measured from recovered map; [5] major 
axis of AMI-LA synthesized beam; [6] minor axis of AMI-LA synthesized 
beam. 



[DCE08]+ r 


2° 
















(arcsec) 


(arcsec) 


202 


3C286 


J0329+2756 


19 


41.4 


28.2 


208 


3C286 


J0329+2756 


16 


52.6 


26.6 


210 ] 

211 J 


> 3C48 


J0329+2756 


19 


48.2 


26.3 


212 ] 












213 


1 3C286 


J0329+2756 


20 


40.2 


29.8 


214 J 












215 


3C48 


J0329+2756 


22 


48.5 


25.9 



Numbers in bold indicate field designations as shown in Fig. 1 



performed using short observations of 3C286, 3C48 and 3C147. 
We assume I-hQ flux densities for these sources in the AMI-LA 
channels consistent with the updated VLA calibration scale (Rick 
Perley, private comm.), see Table 2. Since the AMI-LA measures 
I-i-Q, these flux densities include corrections for the polarization 
of the calibrator sources. A correction is also made for the chang- 
ing air mass over the observation. From other measurements, we 
find the flux calibration is accurate to better than 5 per cent (AMI 
Consortium: Scaife et al. 2008; AMI Consortium: Hurley-Walker 
et al. 2009). Additional phase (secondary) calibration is done us- 
ing interleaved observations of calibrators selected from the Jodrell 
Bank VLA Survey (JVAS; Patnaik et al. 1992). After calibration, 
the phase is generally stable to 5° for channels 4-7, and 10° for 
chaimels 3 and 8. The FWHM of the primary beam of the AMI-LA 
is Ri 6 arcmin at 16 GHz. Due to their superior phase stability only 
channels 4 — 7 were used for this work, resulting in an effective 
total bandwidth of 3 GHz. 

Reduced data were imaged using the AIPS data package. 
CLEAN deconvolution was performed using the task IMAGR which 
applies a differential primary beam correction to the individual fre- 
quency channels to produce the coinbined frequency image. Due 
to the low declination of this sample, uniform visibility weighting 
was used to improve the AMI-LA PSF. The AMI-LA is sensitive 
to angular scales from 0.5 — 6 arcmin, although this varies as a 
function of declination, hour angle coverage and data flagging. 

In what follows we use the convention: 5v °= v", where Sy is 
flux density (rather than flux, Fy = vSv), V is frequency and a is 
the spectral index. All errors quoted are I a. 

Where a known object is not detected by the AMI-LA an up- 
per limit on the 1.8 cm radio flux density of that source is given. 
These linuts take into account the position of the source within 
the AMI-LA primary beam. Since source detection is performed 
in the maps before correction for the primary beam this means 
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Table 1. The AMI-LA sample of embedded protosteUar sources selected from the catalogue of Dunham et al. (2008). Columns are [1] source number from the 
Dunham et al. (2008) catalogue; [2] Right Ascension of source in J2000 coordinates; [3] Declination of soiffce in J2000 coordinates; [4] candidate Group from 
Dunham et al. (2008); [5] IR luminosity from Dunham et al. (2008); [6] identifier from Evans et al. (2009); [7] bolometric luminosity from Evans et al. (2009); 
[8] bolometric temperature from Evans et al. (2009); [9] identification from Harvey et al. (2007); [10] identifier from Enoch et al. (2009); [11] classification 
from Enoch et al. (2009); [12] envelope mass from Enoch et al. (2009); [13] predicted 1.8 cm flux density based on derived correlation between bolometric 
luminosity and radio luminosity (Paper II). 



[DCE08] 


RA 
(J2000) 


Dec. 
(J2000) 


Group 


(L©) 


[EDJ09] 


(Ls) 


Tbol 
(K) 


[HMH07] 


Bolo 


Class 


Menv 
(Mo) 


r-pred 
''l.Scm 

(mJy) 


202 


18 28 44.78 


+00 51 25.9 


3 


0.061 


527 


0.078 


180 


24 


3 


I 


0.24 ±0.02 


0.069 


208 


18 29 09.05 


+00 31 27.8 


1 


0.009 


577 


0.024 


410 


75 


15 





1.16 ±0.02 


0.040 


210 


18 29 49.63 


+01 15 22.0 


2 


0.347 


637 


2.1 


97 


141 


23 





7.98 ±0.07 


0.317 


211 


18 29 51.98 


+01 15 38.2 


3 


0.036 




4.0 


110 






1" 




0.426 


212 


18 29 52.51 


+00 36 11.9 


3 


0.015 


650 


0.85 


100 


154 


24 





0.54 ±0.03 


0.209 


213 


18 29 53.04 


+00 36 06.8 


2 


0.185 


653 


0.57 


860 


157 




1" 




0.173 


214 


18 29 54.31 


+00 36 01.4 


2 


0.105 


661 


0.64 


110 


166 


24 


I 


0.48 ±0.05 


0.182 


215 


18 29 57.67 


+01 13 04.4 


2 


0.018 


676 


0.55 


820 


181 


25/28 


I 




0.170 



" Classification based on the criteria of HatcheU et al. (2007). 



Table 4. AMI-LA observations of the Serpens Core (Cloud A). Columns are [1] AMI-LA source number; [2] Right Ascension of radio source peak; [3] 
Declination of radio source peak; [4] Integrated 1.8 cm flux density; [5] VLA association; [6] Integrated 3.5 cm flux density; [7] Spectral index from 3.5 to 
1.8 cm; [8] Associations from HMH07; [9] Bolometric luminosity from Enoch et al. (2009a), values in itaUcs from Evans et al. (2009); [10] envelope mass 
from Enoch et al. (2009a); & [11] protosteUar classification 



AMI 


RA 


Dec 


'^l.Scmiiiit 


VLA 


'^3.5cm,iiit 


"3.5 


[HMH07] 


^bol 


A'env 


CI. 




(J2000) 


(J2000) 


(mJy) 




(mJy) 






(Lq) 


(Mq) 




DCE08-210: 




















1 


18 29 49.63 


+01 15 21.9 


4.736 + 0.237 


7 


7.54 


-0.73±0.16 


141 


17.3+ 


7.98 ±0.07 







18 29 51.06 


+01 15 33.9 


1.301 ±0.067 


10 


0.14 


0.81 ±0.10 


150 


2.6 


1.23±0.12 


I 


2 






18 29 52.19 


+01 15 47.8 


8 


0.64 












18 29 52.85 


+01 14 56.0 


< 0.101 


11 


0.09 


<0.18 


155 


1.2 


0.74 ±0.07 


I 




18 29 49.13 


+01 16 19.8 


< 0.102 








137 


4.4 


1.53±0.15 


I 




18 29 51.14 


+01 16 40.6 


< 0.110 


9 


0.15 


< -0.49 


146 


1.7 


1.84±0.18 


I 


3 


18 29 48.10 


+01 16 44.9 


0.355 ±0.026 


5 


0.28 


0.38 ±0.44 


135 


2.1 


3.72 ±0.37 





4 


18 29 49.60 


+01 17 07.0 


0.245 ±0.023 








139 


1.1 




I* 




18 29 56.87 


+01 14 46.5 


< 0.125 








176 


5.8 


0.69 ±0.07 


I 




18 29 57.72 


+01 14 05.7 


< 0.148 


16 


0.17 


< -0.22 


182 


13.8 


1.90±0.19 


I 




18 29 58.77 


+01 14 26.2 


< 0.151 








190 


0.62 


2.00 ±0.20 


I 


5" 


18 29 34.63 


+01 15 07.0 


3.139±0.158 


2 


2.17 


0.58 ±0.05 










6" 


18 29 44.05 


+01 19 22.4 


2. 127 ±0.108 


3 


2.10 


0.02 ±0.005 










DCE08-2I5: 




















1 


18 29 56.33 


+01 13 19.4 


0.450 ±0.031 


13 


0.12 


2.09 ±0.00 


175 


2.6 


2.35 ±0.24 


I 




18 29 57.67 


+01 13 04.4 










181 


2.0 


2.56 ±0.03 


I 




18 29 57.89 


+01 12 46.2 




17' 


0.18 




187 


32.0 






8 


18 29 57.84 


+01 12 37.8 


1.000 ±0.050 


17 


1.82 


-1.05 ±0.02 


186 


6.1 








18 29 59.92 


+01 13 11.6 




19 


0.11 




198 


4.0 








18 30 00.70 


+01 13 01.4 




15 


0.10 




203 


1.2 


2.56 ±0.03 







18 29 57.72 


+01 14 05.7 


<0.118 


16 


0.17 




182 


13.8 


1.90±0.19 






18 29 59.56 


+01 11 59.0 


< 0.122 








197 


1.8 


1.14±0.11 






18 29 58.77 


+01 14 26.2 


< 0.128 








190 


0.62 


2.00 ±0.20 






18 30 02.73 


+01 12 28.2 


< 0.128 








208 


2.2 


1.32±0.13 






18 29 56.87 


+01 14 46.5 


< 0.138 








176 


5.8 


0.69 ±0.07 






18 29 52.85 


+01 14 56.0 


< 0.162 


11 


0.09 




155 


1.2 


0.74 ±0.07 





" outside AMI-LA primary beam FWHM. 
* classification from Kaas et al. (2004) 
t Enoch et al. 2009b. 



that the true limit on the unattenuated flux density will be given 
by ^lini < 5 C7th/A(r), where A(r) is the primary beam attenuation 
at a radial distance r from the pointing centre and <Tth is the r.m.s. 
noise on the map. 



4 RESULTS 

In this section we address each of the regions of the Serpens molec- 
ular cloud separately and examine the sources detected in the radio 
data individually. These examinations are made with reference to 
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Table 5. AMI-LA observations of the Serpens Cloud B. Columns are [1] AMI-LA source number; [2] Right Ascension of radio source peak; [3] Declination 
of radio source peak; [4] Integrated 1.8 cm flux density; [5] VLA association; [6] Integrated 3.5 cm flux density; [7] Spectral index from 3.5 to 1.8 cm; [8] 
Associations from HMH07; [9] Bolometric luminosity from Enoch et al. (2009a), values in italics from Evans et al. (2009); [10] envelope mass from Enoch 
et al. (2009a); & [11] protosteUar classification 



AMI RA 


Dec 


•^1.80111,011 


VLA 


'^3.5cm,mt 


"3.5 


[HMH07] 


i-bol 


Mem 


CI. 


(J2000) 


(J2000) 


(mjy) 




(mJy) 






(L0) 


(Mq) 




DCE08-208: 




















9 18 29 00.8 


+00 32 19 


0.141 ±0.023 
















10 18 29 01.7 


+00 29 48 


3.583 ±0.062 


4 


3.30 


0.13 ±0.04 


59 


4.2° 




n/m 


11 18 29 06.4 


+00 30 38 


1.204 ±0.121 


6 


0.32 


2.07 ±0.40 


68 


6.5« 


2.77 ±0.14 





12 18 29 09.1 


+00 31 33 


0.181 ±0.024 








75-A,B,C 


1.7 


1.16 ±0.02 





18 29 02.1 


+00 31 21 


< 0.0703 








60 


0.09 


0.17 ±0.02 


I 


18 29 09.1 


+00 31 33 


< 0.0904 








63 


0.12 


0.22 ±0.02 


I 


DCE08-213: 




















13 18 29 52.2 


+00 38 51.5 


0.589 ±0.067 
















14 18 29 53.3 


+00 34 09.4 


0.267 ±0.082 
















18 29 53.04 


+00 36 06.8 


< 0.0950 








157 


0.57 






18 29 52.52 


+00 36 11.7 


< 0.0950 








154 


1.0 


0.54 ±0.03 





18 29 54.30 


+00 36 01.3 


< 0.0960 








166 


0.17 


0.48 ±0.05 


I 



" from Harvey & Dunham (2009). 



Table 6. AMI-LA observations of the Serpens Filament. Columns are [1] AMI-LA source number; [2] Right Ascension of radio source peak; [3] Declination 
of radio source peak; [4] Integrated 1.8 cm flux density; [5] VLA association; [6] Integrated 3.5 cm flux density; [7] Spectral index from 3.5 to 1.8 cm; [8] 
Associations from HMH07; [9] Bolometric luminosity from Enoch et al. (2009a), values in italics from Evans et al. (2009); [10] envelope mass from Enoch 
et al. (2009a); & [11] protosteUar classification 
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18 28 44.8 


+00 51 26 


< 0.095 




24 


0.04 


0.24 ±0.02 
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18 28 45.0 


+00 52 04 


< 0.098 




26 


1.09 


0.29 ±0.03 


1 


18 28 44.0 


+00 53 38 


< 0.138 




23 


0.19 


0.26 ±0.03 


I 



several archival data sets. Notably these include the Northern VLA 
Sky Survey (NVSS: Condon et al. 1998) at 1.4 GHz in the radio 
regime, the sub-mm catalogue of Enoch et al. (2009) from the Bolo- 
cam instrument at 1.1 mm and the Mid Infra-Red Spitzer catalogue 
of Evans et al. (2009), as well as the Serpens specific Spitzer cata- 
logue of Harvey et al. (2007). 

Source type is characterized by radio spectral index with 
extra-galactic sources expected to have steeply falling spectra with 
frequency, a > —0.5, and thermal sources expected to have shal- 
low or rising spectra —0. 1 ^ a ^ 2. Thermal sources are assumed 
to represent protosteUar objects, where a = —0.1 indicates opti- 
cally thin free-free radio emission and a = 2.0 indicates optically 
thick emission. Values intermediate to these limits are assumed to 
be sources with partially optically thick media (see e.g. Reynolds 
1986; Panagia & Felli 1975). 

These details are consistent with data presented in Pa- 
pers I & II where correlations were seen between the radio luminos- 
ity from objects identified as protosteUar with a number of physical 
characteristics determined from other wavelength regimes. Partic- 
ularly relevant to this work are the correlations with envelope mass 
and bolometric lirarinosity (Paper II) which allow us to make pre- 
dictions of radio limiinosity based on emission in other regimes, 

log[Li.8cm(mJykpc2)] = -(1.74±0.18) 

+(O.51±O.26)log[Lboi(L0)], (1) 



and 

log[Li.8cm(mJykpc2)] = -(2.23 ±0.65) 

+(0.68 ±0.62) log [Menv(MQ)]. (2) 

The predicted radio flux densities for the target sample presented 
here, based on bolometric luminosity, are listed in Table 1. 



4.1 Serpens Core (Cluster A) Region 

Candidate cores DCE08-210, DCE08-211 and DCE08-215 are 
found in the region of the Serpens core. DCE08-210 is associated 
with the sub-mm core Serpens SMM 1 and is closely separated 
from DCE08-211. At high resolution Serpens SMM 1 can be re- 
solved into a triple radio source, with components coincident with 
its core and outflows (Rodriguez et al. 1989; Curiel et al. 1993) 
and is the most luminous object in this molecular cloud. DCE08- 
211 is associated with Serpens SMM 10, which is an order of 
magnitude smaller in mass than its neighbour. DCE08-215 is lo- 
cated to the south-east of the previous two objects and is associated 
wifli Serpens SMM 4, also known as EC95 (Eiroa & Casali 1992) 
and thought to be the X-ray source Ser-X3 (Smith, Giidel & Benz 
1999). 
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4.1.1 Existing radio observations of cluster A 

The Serpens Cluster A region was previously surveyed at 3.5 cm 
with the VLA by Eiroa et al. (1993). Following the Spitzer survey 
of this region it is now possible to match each of the 3.5 cm de- 
tections with individual protostellar sources. These identifications 
are listed in Table 4. The flux densities have been adjusted to bring 
them on to the same flux scale as the AMI-LA, which follows the 
updated VLA calibration (Rick Perley private comm.). We identify 
radio detections at 8.5 GHz from Eiroa et al. (1993) with protostel- 
lar sources from Harvey et al. (2007) when they are coincident to 
within a radius of 10". The protostellar objects detected at 3.5 cm 
are in general deeply embedded objects, with no coincidence with 
known Class U/IU sources. Within the VLA fields there is only 
one source classified as embedded by Harvey et al. (2007b) that is 
not detected in this sample, HMH07-195. Of the detected sources, 
two are not classified as embedded because they possess a 2MASS 
counterpart: HMH07-146 & -182 but are identified as Class I by 
Enoch et al. (2009); two are unclassified, most likely because both 
lack detections at 70 ;Um: HMH07-182 & -186; HMH07-198 is 
classified as cold and embedded by Harvey et al. (2007b) however, 
it is not identified by Enoch et al. (2009) as being part of a sub-mm 
core. Those sources which are not identified in the sub-mm are not 
necessarily more evolved objects but could instead be examples of 
deeply embedded sources with envelope masses below the detec- 
tion threshold of current sub-mm surveys. 



4.1.2 DCE08-210 

The radio source detected here towards DCE08-210 is coincident 
with IRAS 18273+0113, also known as Serpens FIRSl and Ser- 
pens SMM 1. In the mid-infrared this object can be seen to en- 
compass at least five separate sources. DCE08-210 is coincident 
with the near-infared source EC4I (Eiroa & Casali 1992), which is 
also the ISOCAM source K258a (Kaas et al. 2004). This object is 
thought to be a Class I source based on its 0!|^|^''* index and this 
classification agrees with the bolometric luminosity and tempera- 
ture found by Spitzer and given in Evans et al. (2009), where it 
is designated [EJD2009]-637. However, using combined data from 
2MASS, Spitzer and Bolocam, Enoch et al. (2009; 2009b) clas- 
sify FIRSl as a Class YSO with a bolometric temperature of 
56± I2K, a bolometric luminosity of II ±6L0 and an envelope 
mass of 7.98 ± 0.07 Mq, consistent with its original IRAS classi- 
fication (Hurt & Barsony 1996). This is also supported by its ob- 
served molecular outflows, which are highly extended spatially and 
energetic (Davis et al. 1999; Graves et al. 2010), both properties 
more normally associated with Class objects. 

The closely associated object K258b, 13" to the south-west 
of K258a, and detected in the survey of Serpens by Harvey et al. 
(2007) in four Spitzer bands has been suggested to be simply a knot 
of extended emission (Kaas et al. 2004). 

The dust temperatures found from far-infrared and sub-irmi 
studies of DCE08-2I0 (Hurt & Barsony 1996; Davis et al. 1999) are 
30 — 40 K towards FIRSl. This is high for Class objects, which 
typically have 10 < Td < 20K (see e.g. Hatchell et al. 2007; Visser 
et al. 2002). It is uncertain therefore whether Serpens FIRSl, and 
therefore DCE08-2I0, is a late stage Class YSO or an early stage 
Class I YSO. 

DCE08-21I (Serpens SMM 10) lies to the north-east of 
DCE08-2I0, see Fig. 2. At such close separation the 16 GHz data 
is still heavily dominated by the strong emission from DCE08-210 
and it is not possible to attribute the eastern extension of this emis- 



sion conclusively to DCE08-21I rather than the source VLA-8, 
which is positioned intermediate to DCE08-210 and DCE08-21I. 

To the north of DCE08-210 a spur of radio emission is ob- 
served at 16 GHz. There are two distinct peaks within this spur 
which we attribute to HMH07-135 (Serpens SMM 9; VLA-5) 
and HMH07-139. HMH07-135 was identified as Class I by Kaas 
et al (2004) and later as Class by Enoch et al. (2009). HMH07- 
139 (K254) was also identified as Class I by Kaas et al. (2004). 
HMH07-146 (VLA-9), found slightly to the east of these objects, is 
associated with Serpens SMM 5; this source is tmdetected at 1.8 cm 
with the AMI-LA. 

At the eastern edge of this field there is a 16 GHz radio point 
source which is not identified as protostellar in HMH07 or Evans 
et al. (2009). This object may be extra-galactic, however it is coin- 
cident with Serpens SMM 8, which is also coincident with a peak in 
CO emission (Davis et al. 1999) and could therefore be a possible 
radio protostar. 

To the west of the AMI field is the bright radio source 
NVSS J 182934-1-0 11504. There is some debate over the nature of 
this object. Although it is coincident with a 2MASS point source 
it has no Spitzer counterpart. Eiroa et al. (1993) suggest that it is 
in fact extra-galactic, and indeed the flat spectral index measured 
between 3.5 and 1.8 cm may be misleading for this object as the 
1 .4 GHz flux density is dependent on resolution, varying by a fac- 
tor of 1.5 between the NVSS catalogue (Condon et al. 1998) and 
the higher resolution observations of Snell & Bally (1986). This 
could be evidence of variability, but is more likely to arise as a re- 
sult of flux lost from extended structure at high resolution as the 
object does appear extended in the AMI-LA 16 GHz map. From 
the AMI-LA and NVSS data, which have more comparable visibil- 
ity coverage than the AMI-LA and Eiroa et al. (1993) data, we find 
a[*4 = 0.58 ±0.03. 

4.1.3 DCE08-215 

The region around DCE08-215 contains several radio sources and 
protostellar sources, see Fig. 3(a). We observe a large kidney 
shaped extended source coincident with DCE08-215. This source 
is elongated north-south; Eiroa et al. (1993) detect a number of re- 
solved sources at higher resolution at 8.35 GHz within this area and 
it is likely that the 16 GHz source includes multiple contributions. 
The peak of the 16 GHz radio emission is coincident with VLA- 
17 (Eiroa et al. 1993) and we can assimie that this source provides 
the dominant contribution to the emission. In order to examine the 
area in more detail we subtract a point source at this position di- 
rectly from the visibility data with a flux density consistent with 
the 16 GHz peak value. After removing this source the remaining 
emission has the morphology shown in Fig. 3(b). A distinct point 
source is seen to the north-west of the pointing centre, coincident 
with Serpens SMM 4 (VLA-13); to the south-east an extended re- 
gion of emission is present, coinciding with the sub-mm emission 
joining Serpens SMM 2 and SMM II. No 16 GHz radio emission 
is seen towards Serpens SMM 3. 

We associate the northern point source with K308, which was 
identified by Kaas et al. (2004) as a Class I object. It is not in- 
cluded in HMH07 due to the lack of a detection at 24 flm. Con- 
versely, Winston et al. (2007) identify this object as Class 11, how- 
ever here we retain the Class I designation due to the presence of a 
sub-mm core and we identify this source with Ser-emb 22 (Enoch 
et al. 2009). The extended emission to the south of DCE08-215 
is more difficult to divide conclusively. We tentatively divide this 
emission between the two Class I sources K33I (Kaas et al. 2004) 
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(a) (b) 

Figure 3. (a) DCE08-215 field; (b) DCE08-215 field with point source subtraction. Grayscale and contours are from the AMI-LA at 16 GHz, with contours at 
5,6,7,8 <T etc where a = 22 ;i Jy beam^ ' . Crosses {'+') mark the positions of Class protostellar objects, stars ('*') mark the positions of Class I objects, see 
Table 4. The position of candidate embedded objects from DCE08 are shown as squares. The AMI-LA primary beam FWHM is shown as a solid circle and 
the synthesized beam as a filled ellipse in the bottom left comer. The area shown in (b) is marked as a rectangle in (a). Thick contours in (b) are Scuba 850 ^im 
data from Di Francesco et al. (2008), shown in increments of 2 per cent from a level of 15 per cent. 




1B'>30'"00= 29"50' a9°40= 

Right Ascension (J2000) 



Figure 2. DCE08-210. Grey scale and contours are from the AMI-LA at 
16GHz, with contours at 5,6,7,8 CT etc where fj = 19^iJybeam"'. Crosses 
('-I-') mark the positions of Class protostellar objects, stars ('*') mark the 
positions of Class I objects, see Table 4. The position of candidate embed- 
ded objects from DCE08 are shown as squares. The AMI-LA primary beam 
FWHM is shown as a solid circle and the synthesized beam as a filled ellipse 
in the bottom left comer 



associated with VLA-19 (Eiroa et al. 1993) and HMH07-198, and 
K330 (Kaas et al. 2004) associated with HMH07-197. 

4.2 Cluster B Region 

4.2.1 Existing radio observations of Serpens Cluster B 

Observations of the Serpens Cluster B region were made with the 
VLA at 3.5 cm by Djupvik et al. (2006). Of nine radio detections we 
associate four with protostellar objects from Harvey et al. (2007). 
As with Cluster A the majority (three) of these radio sources cor- 
respond to Class O/I objects; the remaining detection is associated 
with HMH07-59, which is identified as Class II/III (Harvey et al. 
2007; Harvey & Dunham 2009). 

4.2.2 DCE08-208 

The object DCE08-208 was identified as an external galaxy can- 
didate by Evans et al. (2009; [EDJ2009]-577) with a bolometric 
luminosity of 0.024 Lq, but is considered a Group 1 candidate for 
an embedded protostellar source by DCE08. Harvey et al. (2007a) 
who initially identify this object as HMH07-75, divide it further 
into three components A, B and C (Harvey et al. 2007b). Of these 
components they make preliminary classifications based on bolo- 
metric temperature of Class I, I and 0, respectively. Component 
C has a 70 flm counterpart, whereas components A and B do not, 
and it has been suggested that, although component A is point-like, 
component B may simply be extended emission from the outflow 
of component C. The combined bolometric luminosity of the three 
sources is heavily dominated by the Class component, C, which 
has a re-derived value of 1.6 L© (Harvey & Dunham 2009). 
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Radio emission from DCE08-208 is detected at 16 GHz as an 
extension to a brighter adjacent source, see Fig 4(a). It is likely 
that this emission is dominated by HMH07-75C, and appears to be 
extended towards the north of this source. There is no correspond- 
ing radio source detected at 3.5 cm (Djupvik et al. 2006). Assum- 
ing that DCE08-208 was not significantly removed from the phase 
centre of the 3.5 cm VLA observations this would give the radio 
emission a spectral index of a ^ 1.13. 

Immediately to the south-west of DCE08-208 is a bright ra- 
dio source coincident with objects HMH07-67 & -68. We associate 
this source with the Class object HMH07-68 (Harvey & Dun- 
ham 2009), based the on proximity of the radio peak to the Spitzer 
position, however it may also have a contribution from the closely 
removed Class I YSO HMH07-67. HMH07-68 was also detected 
at 3.5 cm by Djupvik et al. (2006) and using their flux density we 
measure a spectral index for this object of a = 2.07 ±0.40. Such 
a steep rising spectral index is typical of optically thick free-free 
emission, such as is seen towards compact Hll regions. These types 
of massive HlI regions are not typical for Class protostellar ob- 
jects and it is therefore more likely that the significant increase in 
flux between 3.5 cm and 1.8 cm is due to the presence of either 
extended structure around this object which has been lost in the 
higher resolution VLA observations, or alternatively variability in 
the radio emission from HMH07-068. 

Further to the south-west there is a third radio source, which 
is centred on the two sources HMH07-58 and -59. HMH07-59 
was also detected at 3.5 cm by Djupvik et al. (2006) and it was 
suggested that the emission from this more evolved Class II/III 
protostellar object could be due to gyro-synchrotron emission. If 
we assume that the flux density seen at 16 GHz is entirely due to 
HMH07-59 we find that this object has a slightly rising spectrum 
with a — 0.13 ±0.04. Again such a spectrum could imply free- 
free as the dominant emission mechanism, however it is also not 
inconsistent with the previously proposed emission mechanism of 
gyrosynchrotron and this is discussed further in Section 6.4. The 
elongated nature of the detection suggests that it may also have 
contributions from HMH07-54 & -56, which are again identified 
as Class II/III objects. Fitting the 16 GHz emission as the sum of 
two point sources, the first centred on HMH07-59 and the sec- 
ond on HMH07-56, we can constrain this additional emission to be 
0.22 mjy. The fact that this second component was not detected 
at 3.5 cm implies that it also possesses a rising spectrum. 



4.2.3 DCE08-213 

There are two Class objects close to the pointing centre of the 
DCE08-213 AMI-LA field. Neither of these is detected, and indeed 
no protostellar sources are detected at 16 GHz in this field with flux 
densities > 96 /J Jy. Two radio point sources are detected away from 
the pointing centre. The radio source to the north is identified with 
NVSS J182952-H003857 and has 5i.4ghz = 2.8±0.6mJy. Com- 
bining these data with the AMI-LA 16 GHz flux density shows 
that this is a steep spectrum extra-galactic radio source with a = 
—0.58 ±0.07. The second source has a flux density approximately 
half that of the northern source and, assuming a similar spectral 
index, this would place it below the detection limit in the NVSS 
survey at 1.4 GHz. With no known protostellar object in the vicin- 
ity of this source we identify it as a non-thermal extra-galactic radio 
source. 




le'^gg^oo" 50= 40= 30= 



Right Ascension (J3000) 

Figure 5. DCE08-202 (no protostellar sources detected). Greyscale and 
contours are from the AMI-LA at 16 GHz, with contours at 5,6,7,8 fj etc 
where a = 20;Ujy beam^' . Stars ('*') mark the positions of Class I objects, 
see Table 6. The position of the candidate embedded object from DCE08 is 
shown as a square. No Class objects ai'e present in this field. The AMI-LA 
primary beam FWHM is shown as a solid circle and the synthesized beam 
as a filled ellipse in the bottom left comer 

4.3 Serpens Filament 

4.3.1 DCE08-202 

The DCE08-202 field in the Serpens filament is completely de- 
void of radio sources at 16 GHz, see Fig. 5. The nearest source 
is almost 4 arcmin to the east of the pointing field and is associ- 
ated with NVSS I182901-I-005117. Flux densities measured this 
far outside the AMI-LA primary beam are imprecise due to un- 
certainties in the primary beam model outside the FWHM but this 
source appears to have a flat spectrum 5'i.4ghz = 3.8±0.5mJy, 
■^leCHz = 3.5±0.4mJy assuming 10% calibration uncertainty at 
16 GHz gives a = -0.03 ±0.01. 

The three protostellar objects in this field are all Class I and 
are uniformly not detected. In the case of HMH07-26 this is sur- 
prising as its bolometric luminosity would suggest that it should be 
detectable according to the correlation between radio and bolomet- 
ric luminosity measured in Paper II. However, it is immediately no- 
ticeable from Table 6 that although these objects vary significantly 
in bolometric luminosity they all have similarly low envelope mass 
estimates. In the case of HMH07-26, which when judged on its 
bolometric luminosity would have a 16 GHz flux density of almost 
300 flJy, has only a predicted 16 GHz flux density of 38 ^Jy when 
judged on its envelope mass. 

4.4 Expected contamination by extragalactic radio sources 

At 16 GHz we expect a certain number of extragalactic radio 
sources to be seen within each of our fields. Following Papers I & 
II, to quantify this number we use the 15 GHz source counts model 
from de Zotti et al. (2005) scaled to the IOC survey source counts 
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Figure 4. Seipens Cluster B fields, (a) DCE08-208 field; (b) DCE08-213 field (no protostellar sources detected). Greyscale and contours are from the AMI- 
LA at 16 GHz, with contours at 5,6,7,8 o etc where a = 16/iJy beam"' for (a) and a = 19/Jjy beam"' for (b). Crosses ('+') mark the positions of Class 
protostellar objects, stars ('*') mark the positions of Class I objects, and diamonds {'o') mark the positions of Class 11 sources, see Table 5. The position of 
candidate embedded objects from DCE08 are shown as squares. The AMI-LA primary beam FWHM is shown as a solid circle and the synthesized beam as a 
filled ellipse in the bottom left corner. 



(Davies et aL 2010). The average rms noise from our datasets is 
~ 19/xJybeam"' and from this model we predict that we should 
see 0.07 sources arcmin"^, or ~ 2 radio sources within a 6arcmin 
FWHM primary beam above a. 5 a flux density of 95 fl]y. Within 
the Serpens core region we would therefore expect to detect ~ 4±2 
extragalactic radio sources, and in the Serpens/G3-G6 region a 
further ~ 9 ± 3 sources. Making the assumption that all sources 
which cannot be identified with a previously known protostellar 
object are extragalactic, as are sources with non-thermal spectra, 
we find 3 radio sources in the Serpens Cluster A fields, consis- 
tent with our prediction; in the remaining fields we find 3 fur- 
ther sources, which is lower than predicted, but consistent at 2 CT. 
The low number of extra-galactic sources detected in these fields 
compared to the model may represent an over-estimation from the 
15 GHz source counts, which are extrapolated from a completeness 
limit of 0.5 mjy - high compared with our detection limit of 95 /ijy. 
One break in the observed source counts is already known at mJy 
flux density levels (Davies et al. 2010) and it is quite possible that 
a further break exists below 0.5 mJy. 



4.5 Expected radio flux density from thermal dust emission 

At 16 GHz there is expected to be a small contribution to the ra- 
dio flux density of protostars due to the long wavelength tail of the 
thermal dust emission from the envelopes around these young ob- 
jects, and this contribution may be estimated using sub-mm data. 
However, the complexity and multiplicity of objects in Serpens 
makes accurate flux density calculation from sub-mm maps diffi- 
cult; sources are crowded and several sources with near neighbours 
may be associated with the same sub-mm core. For the purposes 
of this work we therefore constrain the contribution of thermal 



dust emission to the measured 1.8 cm flux densities of our detected 
sources by extrapolating modified greybody spectra for each ob- 
ject directly from 1.1 mm flux densities. As the flux densities them- 
selves are not listed by Enoch et al. (2009a) we recover them from 
the envelope masses in that work for Class and Class 1 objects in 
Serpens, following 



env^ 1 . 1 mm 



1 mm 



£»2 



(3) 



where D = 260 pc is the distance to Serpens, fCi = 
O.OlWcm^g"' is the opacity and 6i.imm(7d) is the Planck func- 
tion at a dust temperature of = 15 K. The values of these pa- 
rameters are chosen to match those used by Enoch et al. (2009a) 
to derive envelope mass estimates from flux densities. The opacity 
is taken from the dust model of Ossenkopf & Henning (1994) for 
grains which have coagulated over 10^ years and which have thin 
ice mantles (OH5; indicating Column [5] of Table 1 in Ossenkopf 
& Henning 1994). 

The source HMH07-139 (AMI-4) has no counterpart in the 
Bolocam catalogue of Enoch et al. (2009a), however it is coinci- 
dent with a distinct sub-mm peak in the higher resolution 1.3 mm 
MAMBO map of Kaas et al. (2004) where it is identified as a 
Class I object. Following the method of Enoch et al. (2009) we ex- 
tract a sub-mm flux density for HMH07-139 by integrating the flux 
density within a 30" aperture, centred on the Spitzer position, from 
the publically available 850 flm data of Di Francesco et al. (2008). 
From these data we measure a flux density, Sgso^m = 1-9 Jy, and 
consequently an envelope mass of Mgnv = 1.16Mr;i, see Table 4. 

Consequently, only one of the Class 0/1 radio detections lacks 
a sub-mm flux density and envelope mass estimate, HMH07-186. 
This source is not identified with a distinct sub-mm peak, but in- 



10 Scaife et al. 



stead is located between several adjacent compact cores. Without 
evidence for a peak at the position of this source we do not attempt 
to fit a flux density and envelope mass. 

From the 1.1mm flux densities we estimate the thermal 
dust contribution at 16 GHz by extrapolating a modified greybody 
model of the form Sy = v'^i?v(rd) from the 1.1 mm measurement. 
To remain consistent with the envelope mass estimates we again use 
a dust temperature of = 15 K and an index of jSoHS = 1-85. Us- 
ing an alternative dust temperature of 10 K would increase the flux 
densities by approximately 30%. The predictions themselves have 
uncertainties of 20-30% from the errors on the 1.1 mm flux densi- 
ties alone. We note also that, as indicated in Ossenkopf & Henning 
(1994), it is unlikely that j3 will remain constant from 1.1 mm to 
1.8 cm. However in the absence of additional data it is not currently 
possible to constrain this model further. 

In what follows we use the measured flux density at 1.8 cm 
with the predicted greybody contribution subtracted to calculate the 
radio luminosity. This is to ensure that the values used are represen- 
tative only of the radio emission and do not include contributions 
from the thermal dust tail which varies greatly between sources and 
which might therefore influence any conclusions being drawn from 
the correlations examined in the later stages of this paper. In gen- 
eral the corrections made to the radio flux densities for this sample 
are small, in the range 1 — 5%. Notable exceptions are AMl-3 and 
AMI-7 with corrections of 20 and 10 per cent, respectively. We as- 
simie that any thermal dust contribution at 3.5 cm is negligible and 
make no correction for this quantity in the flux densities taken from 
the literature. 



5 CORRELATIONS 

We consider both the 3.5 cm and 1.8 cm data for these fields when 
forming correlations with the IR and sub-mm characteristics of the 
detected objects. In doing this we make two adjustments to the mea- 
sured data. Firstly, for the 3.5 cm flux densities we apply a scaling 
factor of 1.49 = (3.5/1.8)°-^, in order to compare them directly 
with the 1.8 cm flux densities. A spectral index of a = 0.6 corre- 
sponds to partially optically thick free-free emission from both a 
spherically symmetric density distribution (Panagia & Felli 1975) 
and a coUimated density outflow distribution (Reynolds 1986) and 
is generally assumed for the radio emission from protostellar ob- 
jects. Secondly we scale the envelope mass values from Hatchell 
et al. (2007) used in Paper II by a factor of 0.33, in order to compare 
them directly with the mass values for objects in Serpens which are 
taken from the Bolocam catalogue of Enoch et al. (2009). This fac- 
tor is the average ratio of envelope masses for the sources from 
the Perseus region studied in Paper II when measured with Bolo- 
cam (Enoch et al. 2009) and with SCUBA (Hatchell et al. 2007). 
This factor accounts for differences in the assumed dust tempera- 
ture, Td, the opacity, Ky, after correcting for the assumed distance 
to Perseus, between these works when calculating mass values. We 
note that the limiting mass for the Enoch et al. (2009) Bolocam 
catalogue is Menv ^0.1 M©, significantly below the lower limit in- 
dicated in Figs. 6 & 7. 

It can be seen. Fig. 6, that the 3.5 cm data are consistent with 
the correlations measured in Papers I & II. There are two notable 
outliers in Fig. 6, DCE08-210 and HMH07-59. Bolometric lumi- 
nosity estimates for DCE08-210 (Serpens SMM 1) vary; here we 
adopt a bolometric luminosity for this source of Lboi = 17-3 
(Enoch et al. 2009b). This value is higher than that determined 
solely from broadband sub-nmi flux densities (Enoch et al. 2009a, 



see Table 4) as it is calculated using the radiative transfer model of 
Dullemond & Dominik (2004) and takes into account inclination 
effects. 

DCE08-210 is also detected in the NVSS survey at 1.4 GHz 
(Condon et al. 1998) with a flux density 5i,4ghz = 4.5±0.5mJy 
implying that it has a rising spectrum from 21 cm to 1.8 cm. The 
flux density measured at 3.5 cm is in excess of both the 21 cm and 
1.8 cm values, preventing a single power law fit. This may imply 
that DCE08-210 was "flaring" at the epoch of the 3.5 cm VLA ob- 
servations, which could explain the deviation of this datum from 
the general correlation. HMH07-59 is identified in Harvey & Dun- 
ham (2009) as a Class Il/ni source and consequently is the only 
object of its type included in this sample. It is not certain therefore 
that such an object should follow the general trend observed for 
Class and Class I objects. 



6 DISCUSSION 

6.1 Detections and non-detections 

We show both the detected and undetected sources in Figs. 7(a) & 
(b). In these plots data from Papers I and II are shown as unfilled 
circles, detections from this work are shown as filled circles, and 
non-detections from this work are shown as stars at the position of 
their 5 a detection limits as listed in Tables 4. 5 & 6. The plot area 
below the average 5 tr detection limit at the field centres for the full 
sample of 95 /xJy is shaded in dark grey. The area corresponding to 
the detection limit between the field centres and the primary beam 
FWHM is shaded in light grey. Using the distribution of the data 
from Papers I & II, 1 ct bounds are shown as dashed lines for the 
measured correlations, which are shown as solid lines. 

As has already been remarked upon, in Fig. 7(a) there are 
two detected sources which lie significantly above the bounds on 
the measured correlation. These are DCE08-210 (Serpens SMM 1) 
and HMH07-59. Both of these objects also have discrepant 3.5 cm 
flux densities, see Fig. 6(a). HMH07-59 is a Class II/III source and, 
as the only source of this kind in the data, it is uncertain whether 
HMH07-59 should foUow the general trend derived for Class O/I 
objects. Since the radio emission mechanism for Class II/III objects 
is unknown, and considered to be different from that of Class O/I it 
seems possible that this object could deviate from the correlation. 
The 3.5 cm flux density for DCE08-210 has already been suggested 
to represent a 'flare' in the emission of this object, a possibility that 
is supported by the increased removal of this source from the corre- 
lation at 3.5 cm compared with 1.8 cm. The bolometric luminosity 
used here is Lboi = 17.3 Lq, however significantly higher luminosi- 
ties have been reported for this object which may explain the slight 
offset of this object from the trend. 

In Fig. 7(a), there are a number of undetected sources which 
he below the detection threshold, but outside the bounds on the 
general trend, in the region L(,oi > 2.0 Lq. There are foiu- vmde- 
tected sources in this region; three of these are located in multiple 
fields and in this case the limit corresponding to the field where the 
source is closest to the phase centre is plotted. Noticeably, however, 
this population does not appear correspondingly discrepant when 
considered in terms of their envelope mass. Fig. 7(b). A possible 
reason for this is the high ratio of bolometric luminosity to enve- 
lope mass measured for these objects. Excluding these four objects 
the mean ratio, r = L|,oi//Wenv, for this sample is r= 1.13. IIMH07- 
182 and -176, which have the highest limiinosities of those sources 
that are undetected, have r = 7.26 and r = 8.41, respectively. The 
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l°gio[ Lboi (Lq) ] log, J M^„, (Mg) ] 

Figure 6. (a) Correlation of bolometric luminosity with 1.8 cm radio luminosity; (b) correlation of envelope mass with 1.8 cm radio luminosity. Unfilled circles 
show data from Papers 1 & II in (a) and from Paper 11 in (b). Filled circles show 3.5 cm data from Eiroa et al. and Djupvik et al. 3.5 cm radio flux densities 
have been scaled by a factor of 1.49, see text for details; masses from SCUBA data (H07; Paper II) have been scaled by a factor of 0.33 to match those derived 
from Bolocam (E09), see text for details. Best fitting correlations from Paper II are shown as solid Unes. 




'°9lo[ Lboi (Lq) ] log,„[ M^„^ (Mq) ] 

Figure 7. (a) Correlation of bolometric luminosity with 1.8 cm radio luminosity; (b) correlation of envelope mass with 1.8 cm radio luminosity, where masses 
from SCUBA data (H07; Paper II) have been scaled by a factor of 0.33 to match those derived from Bolocam (E09), see text for details. Data points from 

Papers I & II are shown as open circles; data points from this work are shown as filled circles for detected sources and stars ('*') at the position at the 5 a upper 
limit for non-detections. Best fitting correlations from Paper 11 are shown as solid lines, with the 1 a bounds on these fits indicated by dashed lines. Radio 
luminosities below the detection limit at the phase centre are shaded in dark grey, and the region shaded in light grey indicates the range of radio luminosities 
which He between the detection limit at the phase centre and the detection Hmit at the half power point of the primary beam. 



remaining two objects in this region, HMH07-137 and -208, have 
r = 2.87 and r = 1.67, respectively. 

Neglecting uncertainties in the assumed distance to Serpens, 
which are commented upon later, the errors on radio luminosity 
in these correlations are dominated by the conservative absolute 
calibration error of the AMI-LA of 5 per cent. Uncertainties on the 
envelope mass have contributions from the assumed opacity index 
and temperature, as well as a 20-30 per cent uncertainty from the 
sub-mm flux density, whereas errors on the bolometric luminosity, 
as stated in Enoch et al. (2009) are approximately 10 per cent. 

The relationship of L(,oi and Mgny is often used to describe 
the evolution of embedded YSOs. The Lboi — Mjnv diagram for 
this sample is illustrated in Fig. 8. We make a simple illustration 
of the evolutionary behaviour of the YSOs in this diagram un- 
der the assumption that the accretion rate is proportional to the 
envelope mass, M^cc °= ^env(f) (Bontemps et al. 1996; Motte & 
Andre 2001). In this scenario collapse is initiated within a finite 
reservoir of material with inass Mg„^,, with accretion proceeding 
at a rate proportional to the mass of the residual envelope as a 
function of time, where Mace = M^ml'^ and T is a characteris- 
tic timescale taken in this instance to be 10^ yrs. Folowing Bon- 
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Figure 8. M^nv vs. Lboi diagram for Serpens. Class protostars are shown 
as filled circles and Class I as open squares. Unclassified sources are shown 
as open circles. The four radio non-detections that are discrepant from the 
Lrad — Lboi trend, but not the — Me„v trend are shown as stars. Evolu- 
tionary tracks are shown for time dependent accretion, « Menv(0' as 
described in the text, for initial mass reservoirs of 0.2 — 3.2 Mq, as indicated 
on the figure. 
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Table 7. Summary of detection statistics for Tables 4, 5 & 6. The original 
sample includes only objects from DCE08, listed in Table 1. The extended 
sample includes all sources from E()9. inchisive of the original sample. 



Class 


Present 


Detected 


% 


Original sample: 






VeLLO 








0/100 





3 


2 


67 


I 


4 


1 


25 


Group 1 


1 


1 


100 


Group 2 


4 


1 


25 


Group 3 


3 








Extended sample: 









6 


4 


67 


I 


17 


3 


18 
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Figure 9. Correlation of radio luminosity with outflow momentum flux. 
Data from Papers 1 & 11 are shown as unfilled circles, data from this worli 
are shown as filled circles. The predicted relationship between momentum 
flux and radio luminosity (Curiel et al. 1989) for a 100% efficiency and 
a temperature of 3000 K (dot-dash line), lO"* K (sohd hne) and lO' K is 
indicated. 



temps et al. (1996), the envelope mass is assumed to be expo- 
nentially decreasing and consequently the accretion rate has the 
functional form Mace = {Meay/'u)exp—t/'Z, and the stellar mass: 
j»^* = (jW^nv)(l-exp-f/T). 

Under a simple standard stellar model (e.g. Adams et al. 1987) 
we may assume that during the infall stage the bolometric luminos- 
ity is heavily dominated by accretion, Lboi = ^acc = GM^Macc/^*- 
We adopt a value of /?* = 3 R© (Stabler 1988), but note that for a 
more sophisticated model the stellar radius will in fact be a func- 
tion of both the accretion rate and the stellar mass, as will the stel- 
lar luminosity, which will also contribute an additive term to the 
bolometric luminosity such that Lbol = ^acc However in this 
simple scenario where we limit ourselves to the earlier stages of 
evolution we consider this term to be negligible. Models which are 
required to provide information past the characteristic time-scale 
will require these additional terms (see e.g. Motte & Andre 2001). 

Compared to the evolutionary diagrams for the Taurus and p- 
Ophiucus molecular clouds (Motte & Andre 2001) the protostellar 
objects in Serpens seem to have evolved from initial mass reser- 
voirs which are a factor of 2 larger. This apparent difference 
might be explained by the the fact that the entire envelope will not 
be completely accreted onto the central star, but that a proportion 
of it will instead be dispersed by the act of outflows (Ladd et al. 
1998). The four objects which are discrepant with the L,-a(j — L|,oi 
trend, but not the L^^^ —M^nv, have values of r which place them 
in a similar region of the evolutionary diagram. Under our simple 
scenario they appear to be objects which have evolved from a rela- 
tively large initial mass reservoir and are approaching their charac- 
teristic timescale, where ^ 50% of their initial mass will have been 
accreted. 

The detection statistics for all objects included in Tables 4, 5 
and 6 which lie within the FWHM of the AMI-LA primary beam 
and have an identified protostellar counterpart are summarized in 
Table 7. The small sky area, and consequently the small number of 
objects of each type present overall, means that no strong statisti- 
cal conclusions can be drawn from this sample with the exception 
perhaps of the low Class I detection rate for the extended sample. 
However, we note that the detection rates for Class and Class I ob- 
jects in the extended sample are similar to those fotmd for Perseus 
(Paper II) of 72 and 27%, respectively. 



6.2 Outflow momentum fluxes 

Graves et al. (2010) calculate outflow momentum values for objects 
associated with four of the sub-mm cores in Serpens. Assuming 
a uniform dynamical age for each of 3 x 10^ years (Davis et al. 
1999) we can estimate the momentum flux for these sources as 
^out = P/t. The highly confused nature of this region precludes 
the use of more sophisticated momentum flux estimation methods. 
We use the maximum measured momentum from the blue and red 
outflow for each object, correcting the values in Table 3 of Graves 
et al. (2010) by a factor of 2 to account for the full velocity range, 
as described in that paper; by a further factor of 10 to account for 
optical depth effects following Bontemps et al. (1996) and Visser 
et al. (2002); and by a small factor to account for the difference in 
distance assumed to Serpens. A comparison of these values with the 
overall correlation of radio luminosity to outflow momentum flux 
from Papers I & n is shown in Fig. 9. As already described in Pa- 
pers I & n, following the model of Curiel et al. (1989) a significant 
proportion of objects exhibit radio luminosities which are incon- 
sistent with free-free emission arising solely as a result of shock 
ionisation in their molecular outflows. The four objects measured 
here (SMM 1, 4, 3/6 and 8) also appear to have radio emission in- 
consistent with this model, most markedly in the case of Serpens 
SMM 1. In Fig. 9 we show the limit set by the shock ionisation 
model of Curiel et al. (1989) on radio luminosity for a canonical 
electron temperature of 10^ K, typical of free-free emission, as a 
solid line. Also indicated is the same model for a temperature of 
3000 K, the electron temperature indicated by Herschel observa- 
tions of shocks (van Kempen et al. 2010), and a temperature of 
10^ K. It is evident that even for very high temperature electron gas 
the shock ionisation model cannot account for the radio emission 
from a significant proportion of these objects. 

However, the basis of a time dependent accretion rate, such as 
that outlined in the simple exponential model of § 6.1, is founded 
on the assumption that there is a relationship between outflow force 
and accretion rate. This assumption arises from the hypothesis that 
the energy required to provide the outflow is released through ac- 
cretion onto the central protostar (Bontemps et al. 1996). If this is 
the case, and the observed correlation between radio luminosity, 
SyLp-, and envelope mass shown in Fig. 7 is a true one, then there 
should be a corresponding relationship between outflow momen- 
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turn and radio luminosity. Therefore it seems likely that the lack 
of a strong correlation in Fig. 9 is due to the inconsistent methods 
used to derive momentum flux across the sample. 

6.3 Distance dependence 

In the derived correlations and physical characteristics we have 
used a distance for the Serpens molecular cloud of 260 pc, consis- 
tent with the Spitzer catalogue. However, recently Dzib et al. (2010) 
have suggested that Serpens has a distance of 415 pc instead, using 
parallax measurements with the VLBA. This difference in assumed 
distance will affect the luminosity, mass and outflow momentum 
flux values quoted here. Radio luminosities will be raised by a fac- 
tor of approximately 2.5, as will bolometric luminosities and en- 
velope mass values, which are proportional to LP-. Since each of 
these quantities is increased proportionately to the other no change 
will occur in the degree of correlation shown in Figs. 6 & 7. Out- 
flow momentum flux is linearly proportional to distance and con- 
sequently, the outflow momentum flux values shown in Fig. 9 will 
be increased by a factor of 1.6, however the radio luminosity val- 
ues will also be increased by a factor of 2.5. Since the increase in 
radio luminosity is proportionately greater than that of momentum 
flux this acts to worsen the agreement with the model (Curiel et al. 
1989) further. 

6.4 Radio spectral indices 

The spectral index of partially optically thick free-free emission 
from fractionally ionized stellar winds and jets has been shown to 
vary depending on morphology (Reynolds 1986). Shallower spec- 
tral indices {a = 0.25) are found for collimated outflows than in the 
spherical case (0.6 < a < 0.8), and steeper spectra still are found 
for conical jets (a > 0.9). In observations such as those presented 
here it is expected that unresolved contributions will be present in 
the radio spectrum from the core itself as well as a confined jet close 
to this core which then becomes unconfined at larger distances. 

The mean value of spectral index for those sources which are 
detected at both 3.5 and 1.8 cm is | = 0.53 ± 1.14, consistent 
with the canonical value for a partially optically thick spherical 
or collimated steUar wind of a = 0.6. However the possibility of 
variable radio emission from these sources may limit the meaning- 
fulness of this agreement. Certainly in the case of Serpens SMM 1 
(DCE08-210), although we observe a non-thermal spectrum from 

3.5 to 1.8 cm, previous observations at a similar epoch measured 
an optically thin free-free spectrum. Across the AMI-LA band we 
see a similar behaviour with Wami = ^0. 15 ± 0.23, consistent with 
optically thin free-free, although the short frequency coverage of 
the AMI-LA makes the uncertainty on this measurement large. 

We do not fit spectral indices for AMl-7 and AMl-8 in the 
DCE08-215 field as the sources cannot be separated without the 
subtraction of AMl-8 from the map. This subtraction in itself will 
affect the spectrum of AMI-7. For the three radio detections in the 
DCE08-208 field we fit spectra across the AMI band, however the 
low signal to noise for these objects, combined with the short fre- 
quency coverage of AMI causes the uncertainties on these measure- 
ments to be large and consequently we suggest that they should 
only be used as an indication of emission mechanism rather than 
for quantitative comparisons. 

HMH07-59 (AMl-10) has a falling spectral index across the 
AMI band, with Wami = — 135 ± 0.34, differing from the opti- 
cally thin index from 3.5 to 1.8 cm. There are a number of pos- 
sible explanations for this difference, notably the 3.5 cm detection 



is associated solely to HMH07-59 whereas, although we also as- 
sociate AMI- 10 with this source, the 1.8 cm emission is extended 
and probably has contributions from multiple sources. The exten- 
sion of the emission itself will cause a steepening of the spectrum 
across the AMI band, as even a small amount of flux loss over such 
a short frequency range can cause significantly steepening. How- 
ever, HMH07-59 is identified as a Class WOl object and there- 
fore it is possible that the emission we are detecting is not thermal 
bremsstrahlung but gyrosynchrotron. Gyrosynchrotron emission is 
characterised by a sharply peaked spectrum, where the frequency 
of the peak is strongly dependent on the magnitude of the magnetic 
field, B. Following Dulk, Melrose & White (1979) the position of 
this peak can be determined by 
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where is also a weak dependency on the electron density, ne, and 
the path length, d, as well as the electron temperature of the plasma 
and the inclination angle between the wave-normal direction and 
the magnetic field, 6. A flat spectral index from 3.5 to 1.8 cm and 
a falUng spectral index across the AMI band implies that the peak 
of such a gjTO-synchrotron spectrum would be found in the region 
of V ~ 12 GHz. The electron density in such objects is assumed to 
be in the range lO"' ^ "e ^ 10^ cm^-', and the path length, roughly 
equivalent to the infall radius, to be in the range 0.01 ^ d ^ 0. 1 pc. 
Making the assumptions of = 10^ cm"^ and d = 0.1 pc, as well 
as assuming an electron temperature of 10^ K and an inclination 
angle of = ;r/4, this would imply a magnetic field strength of 
B = 270 G, with approximately one order of magnitude uncertainty 
from the possible variation in parameters, as well as a small vari- 
ation in the exponents dependent on the exact optical depth (Dulk, 
Melrose & White 1979). 

Such a magnetic field is not unlikely in the vicinity of a Class 
II/III object. Typical magnetic field strengths of 10s of Gauss have 
been suggested at a distance of a few astronomical units for objects 
with circumstellar disks (Shu et al. 2007), and much higher fields of 
a few kiloGauss have been measured from the surfaces of classical 
T-Tauri stars (Donati et al. 2005; Johns-Krull, Valenti & Koresko 
1999; Johns-Krull 2007). More recently kiloGauss fields have also 
been measured from a Class I object (WL 17; Johns-Krull et al. 
2009). A hypothesis of gyrosynchrotron emission from this object 
would be well tested by full radio polarization measurements, as 
gyrosynchrotron emission is known to be highly circularly polar- 
ized. As the AMI-LA measures I-i-Q it is unable to provide any 
information on the Stokes V component towards this object. 

The remaining two sources, HMH07-68 (AMl-10) and 
HMH07-75 (AMI-11), both appear to have optically thin spectra 
with aAMi = -0. 17 ± 1.08 and ttAMi = — 0.05 ± 1.15, respectively. 
In the case of HMH07-68 this may indicate that the optical depth 
of the ionized medium has reached unity between 3.5 and 1.8 cm, 
typical of ultra-compact Hll. 

The conclusions drawn from comparisons of | and aAMI 
are questionable as the 3.5 and 1.8 cm observations have been made 
at widely different epochs. The question of variability in the radio 
emission of protostars is still unanswered and any variation over 
the period intervening these measurements will make these com- 
parisons invalid. 
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Figure 10. Radio flux density of Serpens SMM 1 adjusted to 5 GHz, see 
text for details, at different epochs. 

6.5 Evidence for radio variability in DCE08-210 
(Serpens SMM 1) 

Serpens SMM 1 was previously observed at 15 GHz by Snell & 
Bally (1986) who found a flux density of 5i5 = 10±3mJy. Dis- 
regarding the large uncertainty in this measurement, it differs from 
the AMI-LA 16 GHz flux density by a factor of ^ 2. The 5 GHz flux 
density measured by Rodriguez et al. (1989) is — 6.4 ± 0.4 mjy, 
a factor of 1.5 lower than that of Snell & Bally (1986). Addi- 
tional evidence for variability may also be identified when com- 
paring the flux densities from Snell & Bally (1986) with that at 
8.5 GHz from Eiroa et al. (1993). Although the spectral index, 
aj^ = 0.0 ±0.3, measured by Snell & Bally (1986) from obser- 
vations made at a similar epoch implies optically thin free-free 
emission, the 8.5 GHz flux density is discrepant from this model by 
5 <7. Similarly the 25.4 GHz flux density of Ungerechts & Gusten 
(1984) of 525.4 = 18 ± 2 mJy is significantly removed from the pre- 
dictions of this model, in the opposite sense, by a degree which 
cannot be accounted for by an enhanced thermal dust contribution 
which is < 1 mJy. 

Evidence for variable accretion in DCE08-210 has previously 
been suggested by Enoch et al. (2009b) who detect a high mass 
disk around this source, which should be unstable and therefore ex- 
hibit luminosity variations. The radio data for this object suggest 
that there was previously a radio "flare" event, the additional emis- 
sion from which has now faded. We show this possibility in Fig. 10, 
where we have assumed that the spectrum of DCE08-210 is opti- 
cally thin from 1 .4 to 22 GHz in a single epoch and have corrected 
the measured flux densities to a common frequency of 5 GHz using 
a canonical spectral index for optically thin free-free emission of 
a = -0.1. 



7 CONCLUSIONS 

The aim of this sample was to increase the data for Class I objects 
and improve the mass to radio luminosity correlation measured in 
Paper II. The low detection rate for Class I objects in Perseus is 
however echoed in Serpens and only 18% of Class I objects present 
have been detected (3 sources). Although this is only a small in- 
crease in the available data these objects do not demonstrate the 
same evolutionary divide hinted at in the correlation of radio flux 
density with envelope mass seen in Paper II. The objects detected 
in Serpens do however follow the general trends measured for ra- 
dio luminosity with both bolometric luminosity and envelope mass, 
and the respective detection rates for Class and Class I objects are 



consistent with those found for Perseus in Paper II. We also con- 
clude that envelope mass provides a better radio detection indicator 
than bolometric luminosity with the deviations from the L^.^^ —M^^nv 
correlation being notably fewer than that of the L^^d — Lboi correla- 
tion. 
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